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Abstract 



We study the real time dynamics of a quantum system with potential barrier 
coupled to a heat-bath environment. Employing the path integral approach 
an evolution equation for the time dependent density matrix is derived. The 
time evolution is evaluated explicitly near the barrier top in the temperature 
region where quantum effects become important. It is shown that there exists 
a quasi-stationary state with a constant flux across the potential barrier. This 
state generalizes the Kramers flux solution of the classical Fokker-Planck equa- 
tion to the quantum regime. In the temperature range explored the quantum 
flux state depends only on the parabolic approximation of the anharmonic 
barrier potential near the top. The parameter range within which the solu- 
tion is valid is investigated in detail. In particular, by matching the flux state 
onto the equilibrium state on one side of the barrier we gain a condition on 
the minimal damping strength. For very high temperatures this condition 
reduces to a known result from classical rate theory. Within the specified 
parameter range the decay rate out of a metastable state is calculated from 
the flux solution. The rate is shown to coincide with the result of purely 
thermodynamic methods. The real time approach presented can be extended 
to lower temperatures and smaller damping. 

PACS numbers: 05.40. -Fj, 82.20.Db 
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I. INTRODUCTION 



Barrier penetration phenomena in quantum systems are ubiquitous in physics and chem- 
istry [1]. Since the reaction coordinates describing the transition across the barrier typically 
interact with a large number of microscopic degrees of freedom, a useful theory has to start 
out from a formulation of quantum mechanics which incorporates the effects of a heat bath 
environment. In the classical region of thermally activated barrier crossings dissipation is 
naturally accounted for by the generalized Langevin equation or related methods. Follow- 
ing Kramers [2] transition rates can then be calculated from a nonequilibrium steady state 
solution describing a constant flux across the potential barrier. 

Unfortunately, at present a correspondingly well-founded approach to dissipative barrier 
transmission processes is not available in the quantum regime. Clearly, the description 
of dissipation within the framework of quantum theory has extensively been discussed in 
the literature. For processes involving tunneling Caldeira and Leggett [3] have shown that 
a path integral formulation is particularly suitable. These studies and extensions thereof 
to flnite temperatures by Larkin and Ovchinnikov [4] and by Grabert, Weiss, and Hanggi 
[5] are based on a thermodynamic method for the calculation of transition rates. In this 
approach pioneered by Langer [6] one determines the free energy of an unstable system by 
means of an analytical continuation and extracts the transition rate from the imaginary 
part. This technique has developed through the work of Stone [7] and Coleman and Callan 
[8]. A related approach based on transition state theory has been developed by Miller 
[9] and extended to dissipative systems by PoUak [10]. However, Hanggi and Hontscha 
[11] have shown explicitly that the two approaches are in fact equivalent. Analytical and 
numerical results on imaginary free energy calculations for dissipative metastable systems 
are summarized in an article by Grabert, Olschowski, and Weiss [12]. While the method 
was found to be very successful in explaining experimental data in well-controlled systems 
[13], its range of validity is not exactly known due to the lack of a flrst principle derivation 
of the rate formulas. 

In this and subsequent articles we re-examine the problem of dissipative barrier pene- 
tration on the basis of a dynamical theory. The general framework is provided by a path 
integral description of the time evolution of the density matrix of dissipative quantum sys- 
tems introduced by Feynman and Vernon [14] and extended to a wider class of useful initial 
conditions by Grabert, Schramm, and Ingold [15]. The theory is also presented in a recent 
book by Weiss [16] and a brief outline is given in section II. Thereby we also introduce 
the basic notation. In the remainder of this paper we apply the theory to a system with 
a potential barrier the height of which is large compared to other relevant energy scales. 
This allows for a semiclassical evaluation of the path integrals. In sections III and IV the 
time dependent semiclassical density matrix is determined explicitly in the high temperature 
region where the nonequilibrium density matrix near the barrier top is only affected by the 
harmonic part of the barrier potential. In section V we then show that the time evolution 
of the density matrix allows for a nonequilibrium flux state which is essentially stationary 
within a large time domain and can be shown to be a quantum generalization of the Kramers 
flux solution. These high temperature results are also part of the thesis by one of us [17]. 

In section VI we use the flux solution to derive the decay rate of a metastable state in the 
region of predominantly thermally activated processes. We obtain the well-known formula 
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for quantum corrections to the classical Arrhenius rate [1,12], however, supplemented by 
a definite condition for the range of damping parameters within which this result is valid. 
Thus, the dynamic approach presented here delineates the parameter region where a purely 
thermodynamic rate calculation suffices. On the other hand, the fiux solution cannot only 
be used to extract transition rates. In fact, it has already been applied to nuclear fission. 
There, the number of evaporated neutrons is related to the time the system needs to get from 
the barrier maximum to the scission point which may be derived from the fiux solution [18]. 
In subsequent papers it will be shown that the approach presented here can be extended to 
lower temperatures where the potential barrier is overcome primarily by tunneling processes. 



II. PATH INTEGRAL DESCRIPTION OF DISSIPATIVE QUANTUM SYSTEMS 

In this section we introduce some key quantities and formulas for quantum Brownian 
motion for later use. The notation follows the review by Grabert et al. [15] to which we 
refer for further details. We also introduce the barrier potential considered in subsequent 
sections and provide a formulation in terms of dimensionless quantities. 



A. Model Hamiltonian and barrier potential 

We investigate the dynamics of a system with a relevant degree of freedom that may be 
visualized as the coordinate g of a Brownian particle of mass M moving in a potential field 
V{q). The stochastic motion of the particle is due to its interaction with a heat bath. The 
entire system is described by the model Hamiltonian [19,3] 

H = Ho + Hn + Hon (1) 

where 

p) 

is the Hamiltonian of the undamped particle, 

n=l 2 V"^" / 

describes the reservoir consisting of N harmonic oscillators, and 

N ^ 

HoR = -qYl ^riXn + Yl Ti ^ (4) 

n=l n=l " n 

introduces the bilinear coupling to the coordinate q of the Brownian particle. The last term 
in (4) compensates for the coupling induced renormalization of the potential V{q). The 
number of the bath oscillators N is assumed to be very large. To get a proper heat bath 
causing dissipation we will perform later the limit iV — )• oo with a quasi-continuous spectrum 
of oscillator frequencies. 
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Specifically, we shall consider systems where V{q) has a barrier. In this paper we pri- 
marily investigate the parameter region where the barrier potential can be approximated 
by an inverted harmonic oscillator potential. We shall see that at lower temperatures an- 
harmonicities of the potential field are always essential. In subsequent articles we extend 
the theory to lower temperatures and investigate systems with arbitrary symmetric barrier 
potentials. Assuming that the barrier top is at g = and V{0) = 0, the barrier potential 
has the general form 



v{q) = -\Muy 




(5) 



where the C2k are dimensionless coefficients and qa is a characteristic length indicating a 
typical distance from the barrier top at which the anharmonic part of the potential becomes 
essential. In particular, we assume C4 > so that the barrier potential becomes broader 
than its harmonic approximation at lower energies ^. Also, we restrict ourselves to systems 
where the potential V{q) does not depend on time. 

A natural quantum mechanical length scale for the system near the barrier top is given 

by 

which is the variance of the coordinate in the ground state of a harmonic oscillator with 
oscillation frequency loq. We make the assumption that for coordinates of order qo the 
harmonic approximation for the potential suffices. This means that qa is large compared to 
the width qo and 

e = qo/qa (7) 
is a small dimensionless parameter which will serve as an expansion parameter in the sequel. 



B. Path integral representation and influence functional 

The time evolution of a general initial state Wq of the entire system composed of the 
Brownian particle and the heat bath reads 

W{t) = ex];>{-iHt/h)Woex];>{iHt/h). (8) 

We shall assume that the state Wq is out of thermal equilibrium due to a preparation 
affecting the degrees of freedom of the Brownian particle only. Then 



^In fact, we may put C4 = 1 tliereby fixing tlie lengtli scale Qa- To make the origin of terms in 
subsequent equations more transparent, we shall keep C4 as a parameter, however. 
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(9) 



where the operators (9^, O'^ act in the Hilbert space of the particle, and 

Wp = Z^^eM-(iH) (10) 

is the equilibrium density matrix of the entire system where the partition function Zfi pro- 
vides the normalization. In an initial state of the form (9) the system and the bath are 
correlated. Hence, the customary assumption that the initial density matrix Wq factorizes 
into the density matrix of the particle and the canonical density matrix of the unperturbed 
heat bath is avoided. Some examples of realistic preparations leading to initial conditions of 
the form (9) are specified in [15]. The simplest case is an initial state which is just the equi- 
librium density matrix Wf). Operators Oj^O'^ projecting onto a certain interval in position 
space provide another example. 

Since we are interested in the dynamics of the particle only, the time evolution of the 
reduced density matrix p[t) = tTRW{t) will be considered, where trR is the trace over the 
reservoir. To eliminate the environmental degrees of freedom it is convenient to employ the 
path integral approach [20,21]. In position representation, the equilibrium density matrix 
reads 

W/siq ,x„,q',x'J = J Vq Dx„ exp (^-^S^[q , x„]^ (11) 

where the functional integral is over all paths q{T), x„[t)^ < t < h/S with q{0) = 
q\x„[0) = x^, and q{Tifi) = q^x„[hl3) = x„. The Euclidian action is given by 



5^[g,x„] = S^[q] + Sl[x„] + S^^[q,x„] (12) 



wi 



ith 



N rhfi /I . 2 1 \ 



n=l 



n=l 



SoR[q,x„] = J2 dr [-a„qx„ + q"^- — ^ • (13) 



\ 2m„c<;2 



The position representation of the time evolution operator K{t) = exp[—iHt /h) of the 
entire system is 

K{qf,x„_^,t,q^,x„.) = JVqVxnexp (^^S[q,x„]^ . (14) 

This functional integral sums over all paths g(s), x„(s), < s < t with q[0) = qi, ^^^(O) = a;„;, 
and qit) = qj, x„[t) = x„_^. The action reads 

S[q, Xn] = So[q] + Snlxn] + SoR[q, x„] (15) 
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Wl 



ith 



So[q] = J^ds Qm^^ - V{q)^ 

SoR[q, Xn] = anqxn - q^- — ^ . (16) 

Combining the path integral representations of the three real and imaginary time propagators 
in (8) and (9), one obtains for the functional integral representation of the reduced density 
matrix 

/'(?/,?/,^) = ^ Jdq,dq,dq dq \{q,,q ,q[,q ) 

xJVqVq'Vq exp{^ {So[q] - So[q']) - ls^[q]] F[q, q' , q] (17) 

where the functional integral is over the set of paths q{s) ^ q' (s) ^ q (t) with 

q{0) = q., q'{Q) = ql q{Q) = q' 
q[t) = qf^ q'{t) = q'^, q [fi/S) = q 

and where 

F[q,q\q] =jdx„_^dx„.dx'„. j Vx„Vx'„ Vx„exp{^{SR[x„] + SoR[q,x„] 

-Sr[x'„] - SoR[q',x'„]) - j{Sl[xn] + So^[q,Xn])} (18) 

is the so-called influence functional, a functional integral over all closed paths x„[s)^ x'^[s)^ 
x„[t) of the environment with 

X„{t) = x'^{t) = Xn^, Xn{0) = Xn{h/3) = Xn., Xn{0) = x'^{0) = x'^.. 

Zr normalizes F so that F = 1 for vanishing interaction, i.e. Zr is the partition function of 
the unperturbed bath. The new normalization factor Z in (17) is given by Z = Z^/Zr. The 
position representations of the preparation operators (9^, O'^ in (9) give rise to a preparation 
function 

\{q,q,q\q') = (19) 
i 

describing the deviation of the initial state Wq from the equilibrium state Wfj. 
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C. Dimensionless formulation and spectral density 



Before we proceed it is convenient to introduce a dimensionless formulation. In the sequel 
all coordinates are scaled with respect to the quantum mechanical length scale qo introduced 
in (6). In particular, the scaled barrier potential then reads 



V{q) 



T 



•^2*; 2k-2 2k-2 



k=2 



(20) 



All frequencies are scaled with respect to the barrier frequency loq and all times with respect 
to lOq^ . The dimensionless imaginary time interval is denoted by 



6 = LOoTi fi 



{2l\ 



and the corresponding scaled imaginary time by cr, < cr < ^. Furthermore, we define sum 
and difference coordinates for the real time paths 



and for later purposes also 



X = q — q\ r = [q -\- q')l2 



X = q — q\ r = [q -\- q')/2 



(22) 



(23) 



for the imaginary time path. 

Finally, we introduce the dimensionless spectral density of the bath oscillators 



TT 



N 



8{lO - LOn) 



(24) 



which contains all relevant information on the heat bath. With the help of the spectral den- 
sity sums over the environmental oscillators may be written as integrals which is convenient 
when performing the limit iV — )• 00 with a quasi-continuous spectrum of the environmental 
oscillators. In this limit we get a proper heat bath causing dissipation. 



D. Propagating function, effective action, and damping kernel 

Now, for the harmonic oscillator model of the reservoir, which is equivalent to linear 
dissipation, the functional integrals occurring in the infiuence functional (18) can be evalu- 
ated exactly. The details of the calculation were given previously [15]. As a result one finds 
that the infiuence of the heat bath is described by an effective action containing a nonlocal 
damping kernel. The functional integral (17) gives for the dimensionless reduced density 
matrix 

p{xf,rf,t) = j dxjdr^dg dq J{xf,rf,t,x^,r^,q,q ) \{x^,r^,q,q ) (25) 
where the propagating function 
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J{xf,rf,t,x,,r,,q,q') = Z ^ jVxVrVq exp Qs[a;, r, g ]^ 



(26) 



depending on the action S[a;,r, given below is a 3-fold path integral over all paths 
x{s), r{s), < s < t in real time with 

a;(0) = Xi, r(0) = r,-, xit) = Xf, rit) = rj 

and over all paths g((T),0 < cr < ^ in imaginary time with 5 (0) = q\ q{0) = q. Hence, 
the trajectories contributing to the propagating function are composed of two paths in real 
time and one in imaginary time. An entire path connects rf with Xf but is interrupted since 
Vi ^ q' and Xi ^ q, in general. These intermediate points are connected by the function 
\{xi,ri,q ,q'). Equation (25) determines the time evolution of the density matrix starting 
from the initial state 

(27) 



p{xf, rf,0) = J dq dq X{xf, rf,q,q )pe{q,q ), 

where pg = tTR[W$) in which W$ is the scaled equilibrium density matrix (11) of the entire 
system. 

The effective action in the propagating function (26) is given by [15] 

re [1^2 _ W / 
T^[x,r,q] = i da -q +V{q) + - dcr'k{cr - a') q {cr)q {a 
Jo \_2 2 Jo 

+ / dcr / dsK* (s — ia) q {a)x{s) 
Jo Jo 

+ / ds [xr — V{r + x/2) + V{r — x/2) — ri^{s)x{s)] 



ds'^{s — s') x{s)r{s') / ds'K'{s — s') x{s)x{s') 

2 Jo 



Here, * denotes complex conjugation. The kernel K[s — ia) may be written as 

K{s — ia) = K'{s — ia) + iK"{s — ia). 
The real part of K[s — ia) is given by 

1 00 



K'{s - la) = - g„{s) exp (z//„cr) 



f28) 



(29) 



(30) 



with the Fourier coefficients 

9n{s) 

where the 



did 



to 



COS o^s 



2mr 
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are dimensionless Matsubara frequencies scaled with loq. The imaginary part of K[s — ia) 
reads 

oo 

K"{s-ia) = - ^ z/„(s)exp(zz/„cr) (33) 



e 

n= — oo 



with the Fourier coefficients 



Us) = 2 T — sin(cu3). (34) 



Furthermore, the so-called damping kernel 



7(s) = 2 / ^cos(cus) (35) 

determines the imaginary part of the kernel (29) for purely real times 

For purely imaginary times < cr < ^ the kernel in the Euclidian part of the effective action 
(28) is given by 

k{(T) = -K{-i(t) + 7(0) : 8{(t) : 
2 

= ^ Cn COS (z/„cr) (37) 

n=l 

where we have introduced the periodically repeated delta function 

oo 

: 8{a) := ^ 8{a-ne). (38) 

n= — oo 

The coefficients 

/■°° dto I(lo) v.. 

Cn = 2 



TT CU^ + Z/^ 

= K\l{K\) (39) 

are related to the Laplace transform ^{z) of the damping kernel 7(3). The kernel k{(T) 
satisfies 

/ dak{a) = 0. (40) 
Jo 

We note that all quantities describing the effect of the heat bath can be expressed in terms 
of the spectral density defined in (24). 
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III. SEMICLASSICAL APPROXIMATION AND PARABOLIC BARRIER 



For an anharmonic potential field the 3-fold path integral (26) cannot be evaluated ex- 
actly. In the following we study the semiclassical approximation which will be seen to suffice 
for small e and coordinates near the barrier top. For high temperatures the anharmonic 
terms in the potential (20) give only small corrections and a simple semiclassical approxi- 
mation for a parabolic barrier with Gaussian fluctuations is appropriate. At the beginning 
of this section we specify the equations of motion for the minimal action paths of the effec- 
tive action. In subsection III B the extremal imaginary time path for high temperatures is 
determined, and in subsection III C the extremal real time paths are evaluated using e as an 
expansion parameter. The minimal effective action and the fluctuations around the extremal 
paths are then determined in section IV leading to the semiclassical time dependent density 
matrix for high temperatures. 



A. Minimal action paths 

Variation of the effective action S[a;,r, introduced in (28) with respect to q leads to 
the equation of motion for the minimal action path in imaginary times 

'i - fda' k{a - a')q [a') - ^^^^ = -i tdsK*{s - ia)x{s). (41) 
Jo dq Jo 

The inhomogeneity on the right hand side couples q (cr) to the real time motion. Variation 
of the effective action S[a;, r, with respect to x and r leads to the equations of motion for 
the minimal action paths in real time 

f + -^ ['ds'-f(s - s')r(s') + {V(r + x/2) + V(r - x/2)} = 
ds Jo 2 dr 

i / ds' K'{s — s')x{s') -\- / daK*{s — ia)q (a). 
Jo Jo 

(42) 



anc 



X-— i'ds'-f{s' - s)x{s') + 2^ {V{r + x/2) + V{r - x/2)} = 0. (43) 



d 



S Js 



The above formulae (25)-(28) for the density matrix and the equations of motion (41)- 
(43) for the minimal action paths hold for any potential held [15]. In the following we shall 
consider explicitly the barrier potential introduced in (5) and (20). 



B. Extremal imaginary time path for high temperatures 



Since we seek a solution of (41) in the flnite time interval < cr < ^, it is convenient to 
employ the Fourier series expansion 
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(li'^)=n J2 ?nexp(z//„cr) . (44) 

n= — oo 

This series for q (cr) continues the path outside the interval < cr < ^ as a periodic path 
with period 6. The continuation causes jump singularities in q{(j) and q{(j) at the endpoints 
which must be taken into account when calculating time derivatives of q (cr). In the following 
we shall show self consistently that for endpoints q^q' of the imaginary time path q{(j) and 
endpoints x,-, Xf of the real time path x[s) that are of order 1 or smaller, the q„ are of order 
1 or smaller for high temperatures. Then, the potential (20) may be approximated by 

V{q) = -\q' + 0{e'). (45) 
Inserting (44) into (41) yields the Fourier representation of the equation of motion 

i^n + Cn - 1) ?n = ^^nX - & + lfn[x{s)] + ign[x{s)]. (46) 

The first term on the right hand side, ip„x^ corresponds to a ■.6[cr): singularity of q{(j) 
according to q (0~) — ? (0"*") = x. The coefficient b is related to a ■.6[cr): singularity of q (cr) 



The functionals 



and 



b = 1(0+) -1(0-) = 1(0+) -lie-). (47) 



fn[x{s)]= fdsUs)x{s) (48) 
Jo 



9n[x{s)] = f dsgn{s)x{s) (49) 
Jo 

describe the coupling to the real time motion. The solution of (46) reads 

qn = u„ {iv„ X -b + ifn[x{s)] + ign[x{s)]) (50) 

with the abbreviation 

U„ = + Cn-1) = + k„|7(kn|) - !)• (51) 

The coefficient b must be determined such that the path q (cr) satisfies the boundary condi- 
tions q (0"*") = q \ q = q ■ Due to the discontinuities of the periodically continued path 
at the endpoints, care must be taken in performing the limit cr — )• 0. From (50) and (51) we 
obtain 

1 °° ' X 1 °° ' X 

q{cr) = -- ~ — exp(zz/„cr) + J2 ~ — " 1) exp(zz/„cr) 



n= — oo '* n= — oo 

oo 



^ ^ri{b - ifn[x{s)] - ign[x{s)]) exp(zz/„cr) (52) 
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where the prime denotes the sum over all elements except n = 0. The first sum gives in the 
limit cr — )• 0^ 



lim - 



±- 



(53) 



The second sum in (52) is regular in the limit cr — )• and vanishes. The third sum is also 
regular and one obtains 



I X 

lim q {a) = -bA + - J2 '^n gn[x{s)] T 7; 

cr^O± U _ Z 



(54) 



where 



1 2 ^ 

n=l 



(55) 



and where the relations (/_„(s) = g„{s) and /_„(s) = — /„(s) have been used with gn{-s)^ 
fn{-s) given in (31) and (34). We note that for a harmonic oscillator A is related to the 
variance of position [15]. For a parabolic barrier, however, A has no obvious physical meaning 
since A < for high temperatures. Taking into account the boundary conditions for the 
periodically continued path q (cr), we obtain from (54) 



A 
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J2 ^ngn[x{s)] 



(56) 



with r defined in (23). Finally, using (50) and (56), the Fourier coefficients q„ of the imagi- 
nary time trajectory read 



ft 

zz/„x + z/„[a;(s)] + Z5'„[a;(s)] + — - — ^ Umgm\x{s)\ 



m = — QO 



(57) 



which determine the minimal action path q (cr) as a function of the endpoints for high enough 
temperatures. This result depends on the as yet undetermined real time path x[s). Below 
we shall show that x[s) remains of order 1 for all < s < t. Thus, (57) confirms our 
assumption made at the beginning of this subsection that for endcoordinates at most of 
order 1 the path q (cr) always remains in the vicinity of the top and is also at most of order 
1. 

When the temperature is lowered, i.e. the dimensionless inverse temperature 6 is in- 
creased, |A| becomes smaller and vanishes for the first time at a critical inverse temperature 

Or 



1 2 ^ 

n=l 



0. 



(58) 



For vanishing damping one has A = — |cot(^/2) and therefore 9c = vr, i.e. Tc = huJo/irk-B 
in dimensional units. As a consequence of (58) the Fourier coefficients q„ in (57) diverge 
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for a purely harmonic barrier when 6c is approached. This divergence corresponds to the 
problem of caustics for a harmonic oscillator [21,22]. Since the harmonic potential is always 
an approximation, anharmonic terms of the barrier potential must be taken into account for 
temperatures near Tc even for coordinates near the barrier top. Hence, the analysis presented 
so far is limited to temperatures above the critical temperature Tc. In an earlier work [22] 
we have investigated the equilibrium density matrix near the top of a weakly anharmonic 
potential barrier for vanishing damping down to temperatures slightly below Tc. In the 
dynamical case with damping the imaginary time part of the propagating function can be 
calculated in a similiar way [23]. 

C. Real time paths 

Let us first consider the equation of motion (42) for the real time path r(s). The inho- 
mogeneity reads 

t fds'K'{s - s')x{s') + f AaK*{s - ia)q{a) = I f As' R{s,s')x{s') + F{s). (59) 
Jo Jo Jo 

Here, we have inserted the result (52) for q (cr) and made use of (57) to obtain the right 
hand side where 

oo 

R{S,S') = K'{S - s') + - Un[gn{s)g„{s') - fn{s)fn{s')] (60) 



and 



with 



F{s) = |rCi(s) - ixC2{s) - tCi{s) ['ds'Ci{s')x{s') (61) 
A A Jo 



oo 



n= — QO 

oo 



^2(3) = ^ VnUnfn{s). (62) 

n= — oo 

In the following we will show that for coordinates x, f of the imaginary time path q{(T) 
and endpoints x,-, Xf of the real time path x{s) that are at most of order 1, the real time 
path r(s) remains for high temperatures also inside this region provided the endpoints r,-, rj 
are also at most of order 1. Anharmonic terms in (42) can then be neglected and we have 

1 d 

-— {V{r + x/2) + V{r - x/2)} = -r + 0{e^). (63) 

2 dr 

Hence, the solution for r(s) is straightforward. It has been shown elsewhere [15] that for a 
harmonic potential it suffices to solve the equation of motion for r(s) only for the real part. 
With r(s) = r'(s) + ir"{s) we have from (42) 
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r' + — ds'7(s - s')r'(s' 
as Jo 



F' 



(64) 



where F' denotes the real part of F given in (61). We introduce the propagator G+(s) of 
the homogeneous equation with the initial conditions G+(0) = 0, G+(0) = 1 which has the 
Laplace transform 



(65) 



The solution of (64) is then obtained as 



r'(s) 



G+U) , 

ds'G+is - 



G+is)- 
s')F'(s' 



G+it) 
G+{s 



G+it) 



G+{t) 

f'ds'G+it - s') F 
Jo 



(66) 



Hence, for endpoints of order 1 or smaller the real time path r(s) is also of order 1 or smaller 
for < s < t and high temperatures as assumed above. When the inverse temperature 
is increased F'[s) = fCi(s)/A grows and diverges at ^ = due to the vanishing of A. 
Therefore, the solution (66) becomes invalid and anharmonicities in the equation of motion 
(42) commence to be important for temperatures near Tc. 

The equation of motion for x[s) is homogeneous and can be shown to be the backward 
equation of the equation of motion for r(s) for vanishing inhomogeneity [15]. Correspond- 
ingly, the solution of (43) reads 




G+{t 



G+it) 



-G+it) 



(67) 



This confirms that for endpoints x,-, Xf of order 1 or smaller x[s) is at most of order 1 as 
assumed above. Nonlinear terms in the equation of motion (43) are at most of order and 
can be neglected. 



IV. SEMICLASSICAL DENSITY MATRIX 

Having evaluated the minimal action paths, we may determine the density matrix in the 
semiclassical approximation by expanding the functional integral about the minimal action 
paths. In subsection IV A we first calculate the minimal effective action and in subsection 
IV B we determine the contribution of the fiuctuations about the minimal action paths. It 
will be seen that above Tc a simple Gaussian approximation for the path integral over the 
fiuctuations about the real time and imaginary time paths suffices. 



A. Minimal effective action 

Inserting the minimal action path in imaginary time (52) and (57) as well as the mini- 
mal action paths in real time (66) and (67) into the effective action (28), one gains after a 
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straightforward but tedious calculation for the minimal effective action for inverse temper- 
atures below 6r the result 



T^{xf,rf,t,x„r„x,r) = T^e{x,r) + T^t{xf,rf,t,x„r„x,r). 



Here, 



t — - + t — X 
2A 2 



(68) 



(69) 



is the well-known minimal imaginary-time action of a damped inverted harmonic oscillator 
where 







9 



J2 ^" (kn|7(kn|) - K 



(70) 



We note that for a harmonic oscillator corresponds to the variance of the momentum 
while for a barrier there is no obvious physical meaning as it is the case with A. 

Since the time dependence of the minimal action paths in real time is determined essen- 
tially by the dynamics at a parabolic barrier, it is advantageous to introduce the functions 



(71) 
(72) 
(73) 



A{t) = --Q{t)G4t) 

where Q{t) denotes the step function and 

S{t) = AG+{t) + G+{t)C+{t) 

with 

through which Tit{xf,rf,t,Xi,ri,x,f) can be expressed conveniently. We note, that for a 
harmonic oscillator Ait) is the imaginary part and S{t) the real part of the time dependent 
position autocorrelation function [15]. With (71) and (72) one has 

T^t{xf,rf,t,x„r„x,f) = 



+rx, 



A{t) S 
'A{t) ~ 2AA{t) 

S \ 



-\- r Xf 



2 A{t) 



A{ty \ , S SA{t) 



A{t) J ^ A AA{t) 



-0 + 







2A{t)) 

S A 

+ 



IXXf (s{t) - ^,S{t) 



s{ty 



A{i) 4A(t)2 \ A2 
A 



5(t) - 



+ A 



m 
m 

Aitf 



2A(t)2 



Ait) 



'Siff 



A2 



Ait) 



A2 J 



-A (^w - n^w; 



(74) 
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B. Quantum Fluctuations 



With the minimal actions (69) and (74) we have found the leading order term of the 
path integral for the propagating function (26). The path integral now reduces to integrals 
over periodic paths ^(0) = ^(t) = and ^'{0) = C'(^) = in real time and j/(0) = y{6) = 
in imaginary time describing the quantum fluctuations around the minimal action paths. 

Firstly, let us consider the real time fluctuations. The contribution of the quantum 
fluctuations around the extremal real time paths can be evaluated in the simple semiclassical 
approximation taking into account Gaussian fluctuations. In the corresponding second order 
variational operator the anharmonicities of the barrier potential are at most of order and 
can therefore be neglected. The linear coupling term with the imaginary time path leads only 
to a shift of the real time fluctuations. The second order variational operator is then simply 
given by that of an inverted harmonic oscillator. Using the propagator A[t) introduced in 
(65), the contribution of the Gaussian fluctuations around each of the two real time paths 
q[s) and q'{s) reads [15] 



jv^ exp 



dsC{s) iC{s)+ / du-f{s-u)C{u)-C{s) 
4: Jo \ Jo 



87r|A(t)| 



(75) 



For high temperatures for which the harmonic approximation for the potential suffices, 
the contribution of the quantum fluctuations to the imaginary time path integral can also 
be approximated by the simple semiclassical approximation. Anharmonicities in the second 
order variational operator can then be neglected. Expanding the quantum fluctuations into 
a Fourier series 

oo 

yi'^)=n J2 ?/nexp(z//„cr) (76) 

n= — oo 

and taking into account the boundary conditions j/(0) = y{6) = 0, one obtains the well- 
known result [15] 

jVy exp ^ dcry{cr) (-y{cr) + dcr' k{(T - a') y{cr') - y{cr] 

U^l^^n. (77) 



When the inverse temperature is increased, the simple semiclassical approximation (77) 
for the imaginary time path integral diverges for inverse temperatures near 6c (A — )• 0). 
This corresponds to the problem of caustics in a harmonic oscillator potential as already 
discussed in subsection IIIB. In the harmonic approximation the amplitude of the fluctu- 
ation mode with eigenvalue A of the second order variational operator becomes arbitrarily 
large for temperatures near Tc. As a consequence, the Gaussian approximation breaks down 
near Tc and anharmonic terms of the barrier potential must be taken into account for the 
temperature region near and below Tc. For vanishing damping we have investigated this in 
detail in [22]. For nonvanishing damping the corresponding analysis is performed in [23]. 

Combining (69) and (77), the semiclassical equilibrium density matrix for inverse tem- 
peratures below 6c may be written as 
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(78) 



Here, Z is a normalization constant and the minimal imaginary time action Ti${x^ r) is given 
by (69). 

Having evaluated the effective action and the fluctuation integral for small e, we gain 
the time dependent semiclassical density matrix for small inverse temperatures up to inverse 
temperatures somewhat below 6c. We flnd 



) \{x,,r,,x,r) (79) 

with the propagating function 

_ 1 1 1 / °° \ 

J(x„ r„ t, r,, r) -j= -j== -I 



exp (^e{x,r) + ]^t{x^,r^,i, 



(80) 



The actions YjQ{x^f) and Yjti^x j^r j ^t^Xi^Vi^x^f) are given by (69) and (74), respectively. 
Within the semiclassical approximation the above formula (79) gives the time evolution of 
the density matrix near the top of a potential barrier starting from an initial state with a 
deviation from thermal equilibrium described by the preparation function A(a;i, r,-, x, f). We 
note that using the short time expansions of A(t) and S'(t) (see [15]) one flnds that (79) 
indeed reduces to (27) for t — )• 0. 



V. STATIONARY FLUX SOLUTION 



Now, we consider a system in a metastable state which may decay by crossing a potential 
barrier. We imagine that the system starts out from a potential well to the left of the barrier. 
Metastability means that the barrier height Vf, is much larger than other relevant energy 
scales of the system such as k^T and /icuo, where huj^ is the excitation energy in the well of 
the inverted potential. The dynamics of the escape process may therefore be determined via 
a semiclassical calculation. In particular, we are interested here in the stationary flux over 
the barrier corresponding to a time independent escape rate. Naturally, a stationary flux 
does not exist for all times when the system is prepared at t = in the metastable well. For 
short times preparation effects are important leading to a transient behavior. For very large 
times depletion of the state near the well bottom causes a flux decreasing in time. A plateau 
region with a quasi-stationary flux over the barrier exists only for intermediate times. We 
shall see that for high enough temperatures the flux state depends only on local properties 
of the barrier potential near the barrier top. The time evolution of an initial state in the 
metastable potential can then be calculated with the propagating function in (79). This will 
be done in this section. Firstly, in V A we introduce the inital preparation. Then, in V B 
we determine the stationary flux solution for high temperatures from the dynamics of the 
corresponding nonequilibrium state. 
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A. Initial Preparation 



The initial nonequilibrium state at time t = is described by the preparation function 

X{x,, r,, X, f) = 6{x, - x)8{r, - f)Q{-r,) (81) 

so that the initial state is a thermal equilibrium state restricted to the left side of the barrier 
only. Then, according to (79), the dynamics is given by 

p{xf,rf,t) = j (\.xArtJ{xf,rf,t,x,,r,,x,,ri)<d{-ri). (82) 

To determine a stationary flux solution we seek for an approximately time independent 
state for intermediate times. As already mentioned, the real time dynamics near the barrier 
top is essentially given in terms of the functions Ait) and S{t) introduced in (71) and 
(72). Evaluating these functions for times larger than 1/cur one gets to leading order an 
exponential growth according to 

Ait) = — I — expiiORt) (83) 



anc 



Sit) = -- cot(^) 1 — — - expiujRt). (84) 

^ ^ 2 ^ 2 ^ 2cuR + 7(cur) + cur7'(cur) ^ ^ ^ ^ 

These functions describe the unbounded motion at the parabolic barrier. Here, '^/'{z) denotes 
the derivative of ^f{z)^ and cur is the Grote-Hynes frequency [24] given by the positive solution 
of io^ + c<;r7(c<;r) = 1. The corrections to (83) and (84) are exponentially decaying in time 
(see [15] for details). In the sequel we consider large times with cur^ ^ 1 so that the 
asymptotic behavior (83), (84) holds. We note that then terms in the effective action (74) 
like A[t) — A[t)'^ /A[t) are exponentially small. For cur^ ^ 1 and f = r,-, x = Xi (74) simplifles 
to read 

T^t{xf,rf,t,x„r,) = T^t{xf,rf,t,x„r„x„r,) = 

A{t) 1 fA{t) S{t) \ 

^'^'^ + "'"'^ W) ^ "'"'m) " \W) ^ J 

-0+ " 



4A(t)2 \ a: 
2^(^)2 + 2^^/ A{t) 



+ + + (85) 



where we have kept all terms with XijAit) since the relevant values of Xi will be seen to 
become of order Ait). 

Hence, the density matrix (82) for times cur^ ^ 1 may be written as 



p[xf,rf,t)= dx^dr^ J{xf,rf,t,x^,ri) Q{-ri) (86) 
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where 

J {x -f J T j ^ t, Xj, Tj) — J {x -f J T j ^ t, Xj, Tj, Xj, Tj) 
1 



(87) 

in which 

T^{xf,rf,t,x„r,) = T^0{x„r,) + T^t{xf,rf,t,x„r,). (88) 
The two parts of this action are defined in (69) and (85). 



B. Stationary flux solution for high temperatures 

We shall see that for high temperatures anharmonicities of the barrier potential are 
not relevant for the fiux state. Since then the effective action is a bilinear function of the 
coordinates, the integrals in (86) are Gaussian and can be evaluated exactly. 

The extremum of the effective action (88) with respect to Xi and r,- lies at 

x^ = -2A{t)xf + 2zA{tf-^ 

r'^ = zSit)xj + Sit)'^. (89) 
Introducing the shifted coordinates 

,-y* . ,-y* . ^ 

f^ = - r°, (90) 

the action may be written as 

t{xf, rf,t, X,, ri) = t{xf, rf,t, r°) + S(x„ f,, t). (91) 

The first term is in fact independent of t and coincides with the imaginary time action of 
the equilibrium density matrix, i.e. 

t{xf, rf,t, r°) = J^eixf, rf) (92) 

where '^^[xf^rf) was introduced in (69). For the time dependent second term in (91) one 
obtains 



S(i;8, f,, — ^ 



-Sit).. s{ty-A' 

A(t) ' ' 4A(t)2 



(93) 



With (92) and (93) the time dependent density matrix (82) may now be written as 

p{xf, rf, t) = peixf, rf) g{xf, rf, t) (94) 
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where 



pe{xf,rf) 



1 1 



is the equihbrium density matrix (78) for temperatures well above Tc, and 



(95) 



— , . , , , / dxi df-j Q(—ri — r°) exp f — Sfx,-, f,-, 
87r|A(t)| J ^ P ^ " " ' 



(96) 



is a "form factor" describing deviations from equilibrium. 

To evaluate (96) we first perform the Gaussian integral over the Xj-coordinate. Complet- 
ing the square one finds 



l{xf,rf,t) 



lAI 



With the transformation 



df,0( 



A 



-Vi — r- ) exp 



^4 5(t)2- A 



2 « 



(97) 



Sit) .A. X 

ri = ^ z — ib{t) Xf 



(98) 



the form factor reads 



y{xf,rf,t) 



s{ty 



^ ^| s{ty - A^ 



(99) 



iFiom (84) we see that for lorO < tt one has S{t) < 0. Now, since A < 0, the function is 
nonvanishing only if 



z < -r-f. 



(100) 



Furthermore, for large times we have S{t)/ S{t) = cur and S{ty / {S{ty — A^) = 1 apart from 
exponentially small corrections. Thus, for cur^ ^ 1 the form factor becomes independent of 
time and reduces to 



-f ( I 2 

dz exp < — - (z — iAuj-^iXf ) 
47r|A| "^-oo 4A 



gn(x f,rf) = 

This combines with (94) to yield the stationary fiux solution 

1 



pn{xf,rf) = peixf.rf) 



47r|A| 



dz exp 



4IAI 



(102) 
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where 



u = —rf -\- i\A\ijjRXf. (103) 

To illustrate the result (102) the diagonal part gfi(0,g) of the form factor (101) is shown 
in Fig. 1 as a function of the scaled coordinate q for ohmic damping with 7 = 3 and three 
different values of the inverse temperature 6. One sees that with increasing 6 the width of 
the nonequilibrium state becomes smaller due to the fact that energy fluctuations of the 
Brownian particle decrease for lower temperatures. We note that for 7 = 3 the function A 
vanishes at the critical inverse temperature 6c = 5.79 .. . 

Finally, let us discuss the range of validity of the time independent flux solution (102). 
Firstly, in view of (93), the integral over the x^-coordinate in the expression (96) for the 
form factor is dominated by values of Xi of order \/—A or smaller. After performing the 
Xj-integration, the relevant coordinate fi/ S{t) in the remaining integral (97) becomes of or- 
der 1 /\/—A or smaller. For high enough temperatures Xi and fi/S{t) are then of order 1 or 
smaller which is consistent with the harmonic approximation. However, with increasing in- 
verse temperature |A| decreases and the relevant values of fi/ S{t) become larger. Therefore, 
the harmonic approximation of the barrier potential which leads to the Gaussian integrals in 
(96) is valid only for high enough temperatures. We will show in [23] that anharmonicities 
of the potential barrier become essential for temperatures near the critical temperature Tc 
where the flrst caustic in the harmonic potential arises. 

Secondly, as mentioned above, the flux state (102) was found only for cur^ ^ 1. There 
is also an upper bound of time due to the fact that the flux solution was obtained by 
the harmonic approximation of the potential. To estimate this upper bound we consider 
coordinates within the width of the diagonal part of the flux solution (102) which is of order 
^1 A|. Anharmonic terms of the potential (20) are negligible when e^q^ <^ 1, i.e. q^ <€i qa 
where q^ is the dimensional coordinate. Hence, inserting the rj^-dependent term of the shift 

r° in (89) for values of rf of order ^J\A\ one gains the condition e^S'(t)^/|A| <^ 1. According 
to (84) for large times S{t) becomes of order (c<;r/2) cot(c<;R^/2) exp(c<;Rt) . Furthermore, as 
far as orders of magnitudes are concerned (c<;r/2) cot(c<;R^/2) is of same order as |A| for all 
temperatures well above Tc. This can be seen by estimating A from (55) in various limits. 
As a consequence, anharmonic terms in the potential can be neglected only if exp(c<;Rt) <^ 
l/e^J\A\. Combining the two bounds we see that a stationary flux solution exists within the 
plateau region 



1 1 /2tan(cuR6i/2) 

1 < exp cuRt < ^ -J y ^ ' \ 104 

e^|A| e V tUR 

For very high temperatures ^ ^ 1, where |A| 1/^, the upper bound of time reduces to 



exp (cuRt) ^ a/^/c which can be written as exp (cur^) ^ qay (3mLOQ/2 in dimensional units. 
For lower temperatures, where |A| is of order 1, one obtains exp (iORt) <^ ?a/<Zo in dimensional 
units. We note that (104) holds only for inverse temperatures 6 < 6c where |A| is flnite. In 
summary the flux solution (102) is valid for high enough temperatures and times within the 
range (104). 
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The result (104) has a simple physical interpretation. The exact propagating function 
of the nonlinear potential gives the probability amplitude to reach arbitrary coordinates 
(xf^Vf) when starting from coordinates (xi^Vi). Now, consider endpoints (xf^Vf) near the 
barrier top. Clearly, the main contribution to the corresponding probability amplitude then 
comes from paths with energy of order of the potential energy at the barrier top. Other 
paths either have smaller energy and do not reach the barrier top or much larger energy 
and are therefore exponentially suppressed. However, a trajectory starting at t = far away 
from the barrier top in the anharmonic range of the potential with an energy of order of the 
potential energy at the barrier top needs at least a time of order ln[qa/ Aq^) to reach a region 
of order Aq^ about the barrier top. Now, in dimensional units the flux state has a width of 
order Aq^ = qo^J\A\. Hence, for times sufficiently smaller than ln[qa/ Aq^) = ln(l/ey^jAj) the 
probability to reach endpoints near the barrier top is basically determined by trajectories 
starting also near the barrier top. The upper bound in (104) gives the time where trajectories 
starting far from the barrier region contribute essentially to the nonstationary state. 



VI. MATCHING TO EQUILIBRIUM STATE IN THE WELL AND DECAY RATE 

With (102) we have found an analytic expression for a stationary flux state of a 
metastable system for temperatures above Tc. The function pn^Xf^rf) describes the re- 
duced density matrix of the system for coordinates in the barrier region and for times 
within the plateau region (104). This nonequilibrium state depends on local properties 
of the metastable potential near the barrier top only. On the other hand the metastable 
system is in thermal equilibrium near the well bottom. This means that the flux solution 
must reduce to the thermal equilibrium state for coordinates qj, q'j on the left side of the 
barrier sufficiently far from the barrier top but at distances much smaller than which is 
the typical distance from the barrier top to the well bottom. In this section we verify this 
condition for the high temperature nonequilibrium state (102) and use the result to derive 
the decay rate of the metastable state in the well. 



A. Matching of flux solution to equilibrium state 

For coordinates qf^ q'^ to the left of the barrier the form factor (102) has to approach 1 
as one moves away from the barrier top. Let us investigate the region of coordinates Xf^rj 
where 



47r|A| 



dz exp 



4IAI 



<exp(-l/e^) < 1. 



(105) 



Here, u = u[xf^rf) is given by (103) and the exponent fi is positive. In the semiclassical 
limit /i can be small. Now, in the halfplane rj^ < of the Xfrf-plane the condition (105) 
defines a region which may be delineated by the relation 



1/2 



1-^/ I ~ 



/ i2 A 2 



(106) 
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On the other hand, the equihbrium density matrix (95) is nonvanishing essentially only for 



\xf\<^fl^. (107) 



Now, in order that the flux solution matches to the equilibrium state within the region of 
coordinates where (95) and (96) are both valid, (106) and (107) must overlap in the XfVf- 
plane for sufficiently large rf < but \rf\ <^ . Thereby, one has to take into account 
that due to the definitions of rf and Xf in (22) one has \xf \ < 2\rf \ for qf^q'f < 0. For 
high temperatures this latter relation is always fulfilled by virtue of (107). Hence, from the 
relations (106) and (107) one obtains 



|A|e2-^ < 1 -cu^^ (108) 



where the //-dependence may be disregarded, since we may choose fi <^ 1. 

The fiux solution (102) describes the density matrix of the metastable system in the 
barrier region only when (108) is satisfied. This condition depends on temperature, anhar- 
monicity parameter, and damping. We note that in dimensionless units is the typical 
barrier height with respect to the well bottom. 

So far, we have written all results in terms of dimensionless quantities according to the 
definitions in subsection II C, except for a few formulas where the change to dimensional 
units was mentioned explicitly in the text. In order to facilitate a comparison with earlier 
results we shall return to dimensional units for the remainder of this section. Then , we 
obtain from (108) in dimensional units the condition 



Here, according to (55) and (70) one has in dimensional unit 



s 



1 °° 1 



and 



where loq denotes the oscillation frequency at the barrier top, Vf, the barrier height with 
respect to the well bottom, and the 

.„ = ^ (112) 

are the Matsubara frequencies. The dimensional Grote-Hynes frequency is determined by 
the largest positive root of the equation 

cur + cur7(cur) - cuq = 0. (113) 



23 



iFiom a physical point of view (109) defines the region where the infiuence of the heat 
bath on the escape dynamics is strong enough to destroy coherence on a length scale smaller 
than the scale where anharmonicities becomes important. Only then are nonequilibrium 
effects localized in coordinate space to the barrier region. For very high temperatures, 
i.e. in the classical region UoTifi <^ 1, one has |A| IjiOQUfi and l/hl3. Then, one 
obtains from (109) for Ohmic dissipation with 7(cu) = 7 the well-known Kramers condition 
[2] k^TuolVh <^ 7 where we have used 1 — cu^/cUq ^f/iOo for small damping. When the 
temperature is lowered |A| decreases and the range of damping where the fiux solution (102) 
is valid becomes larger. 

To make the condition (109) more explicit, especially for lower temperatures, one has to 
specify the damping mechanism. As an example we consider a Drude model with 7(t) = 
7c<;d exp(— cud^). Clearly, in the limit cud ^ ^o,7 the Drude model behaves like an Ohmic 
model except for very short times of order I/lob- The Laplace transform of 7(t) then reads 



Since the condition (109) is relevant only for small damping strength, it suffices to evaluate 
(109) in leading order in 7. Then, from (113) one gets 







Expanding A given in (110) up to first order in 7 yields 




'0 




where 



A'(0) 



dA 

57 



7=0 



1 



n= — oo 



E 



00 




^o)^(^D + kn|) 



U„\IjJy) 




Correspondingly, for we gain from (111) 





where 




m 

57 



7=0 





Combining (115)-(119) the condition (109) takes the form 



7 > 




(120) 



loy) 2Vf, tan(c<;o/i/3/2) 1 + 2k tan(c<;o/i/3/2)(c<;D + c<;o)/cud 



where 
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-| CXD II 

Let us discuss the relation (120) in the limit cud ^ ^o,7- For very high temperatures, i.e. in 
the classical region UDfifi <^ 1, the last factor in (120) approaches 1 and (120) reduces to the 
Kramers condition 7 ^ iOo/Vb/S. For lower temperatures where UDfifi ^ 1, but UoTifi < tt, 
the coefficient n becomes large. In leading order we have [25] 

^^ln(c^_ 



TT 



Hence, (120) gives 

tiLol 1 
2Vfctan(cuo/i/3/2) 1 + 2 ln(cuD/i/3) tan(cuo/i/3/2)/7r ' ^ ^ 

Clearly, the range of 7 where the result (102) is valid increases for lower temperatures as 
already stated above. Finally, we note that the conditions (109) and (120) hold only for 
temperatures above Tc where |A| is finite. 



B. Decay rate of metastable state 

Clearly, the fiux solution (102) contains all relevant informations about the nonequili- 
brium state. Specifically, the steady-state decay rate F of the metastable state is given by 
the momentum expectation value in the fiux state at the barrier top. In dimensional units 
one has 

^ = \i^{m) + mph (124) 

where the expectation value (•)£! is calculated with respect to the dimensional stationary 
fiux solution. Here, the normalization constant Z is determined by the matching of the fiux 
solution onto the thermal equilibrium state inside the well. From (124) one has in coordinate 
representation 

(125) 

xj=0 

Here, Jfi is the total probability fiux at the barrier top g = 0. Since the essential contribu- 
tion to the population in the well comes from the region near the well bottom, Z can be 
approximated by the partition function of a damped harmonic oscillator with frequency c<;„ 
at the well bottom, i.e. 

Here, Vh denotes the barrier height with respect to the well bottom in dimensional units. 
Note that the potential was set to at the barrier top. Inserting (102) for rj^ = and (126) 
into (125) gives in dimensional units 
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r 



LO^ CUR 



2tT LOq 



(, 



,n=l 



n 



oo 



+ K\l{\Vn\)+tol 



) 



exp(-/3Vfc). 



(127) 



We recall that the dimensional Grote-Hynes frequency cur is given by the positive solution 
of io^ + c<;r7(c<;r) = lOq. As is apparent from the Arrhenius exponential factor, the rate (127) 
describes thermally activated transitions across the barrier where the prefactor takes into 
account quantum corrections. The above decay rate is well-known from imaginary time path 
integral methods [12] and equivalent approaches [1,11,26]. These methods consider static 
properties of the metastable state only. Here, the expression (127) was derived as a result of 
a dynamical calculation. This avoids the crucial assumption of the imaginary time method 
that the imaginary part of the analytically continued free energy of the unstable system may 
be interpreted as a decay rate. In this sense, we have put the result (127) on firmer grounds. 
Moreover, we have obtain the condition (108) on the damping range where (127) is valid. 



Starting from a path integral representation of the density matrix of a dissipative quan- 
tum system with a potential barrier we have derived an evolution equation for the density 
matrix in the vicinity of the barrier top. This equation was shown to have a quasi-stationary 
nonequilibrium solution with a constant fiux across the potential barrier. In this paper we 
have restricted ourselves to the temperature region where the barrier is crossed primarily 
by thermally activated processes. In this region the fiux solution was shown to be indepen- 
dent of anharmonicities of the barrier potential. This nonequilibrium state generalizes the 
well-known Kramers solution of the classical Fokker-Planck equation to the region where 
quantum corrections are relevant. 

On one side of the barrier the fiux solution approaches an equilibrium state as one moves 
away from the barrier top. We have obtained a condition on the damping strength which 
ensures that the equilibrium state is reached within the range of validity of the harmonic 
approximation for the barrier potential. Again this condition can be shown to be a gener- 
alization of a corresponding condition known from classical rate theory. The fiux solution 
can then be matched to the equilibrium state in the metastable well which yields the proper 
normalization of the fiux across the barrier. In particular, we have used the normalized 
fiux state to determine the decay rate of the metastable state. The result was shown to be 
identical with the well-known rate formula for thermally activated decay in the presence of 
quantum corrections as derived by purely thermodynamic methods. The dynamical theory 
presented here gives in addition a criterion for the range of damping parameters where this 
result is valid. 

A particularly interesting feature of the general approach advanced in this article is the 
fact that it can be extended both to lower temperatures and smaller damping. This will be 
the subject of subsequent work. 



VII. CONCLUSIONS 



26 



ACKNOWLEDGMENTS 

Financial support was provided by the Deutsche Forschungsgemeinschaft (Bonn) through 
SFB237. One of us (G.-L.I.) was also supported additionally through a Heisenberg fellow- 
ship. 



27 



REFERENCES 



[1] P. Hanggi, P. Talkner, and M. Borkovec, Rev. Mod. Phys. 62, 251 (1990). 
[2] H. A. Kramers, Physica 7, 284 (1948). 

[3] A. 0. Caldeira and A. J. Leggett, Phys. Rev. Lett. 46, 211 (1981); A. 0. Caldeira and 

A. J. Leggett, Ann. Phys. (USA) 149, 374 (1983); 153, 445(E) (1984). 
[4] A. L Larkin and Yu. N. Ovchinnikov, Pis'ma Zh. Eksp. Teor. Fiz. 37, 322 (1983) [ 

Sov. Phys.-JETP 37 , 382 (1983)]; A. L Larkin and Yu. N. Ovchinnikov, Zh. Eksp. 

Teor. Fiz. 86, 719 (1984) [ Sov. Phys.-JETP 59, 420 (1984)]. 
[5] H. Grabert, U. Weiss, and P. Hanggi, Phys. Rev. Lett. 52, 2193 (1984); H. Grabert 

and U. Weiss, Phys. Rev. Lett. 53, 1787 (1984). 
[6] J. S. Langer, Ann. Phys. (N.Y.) 41, 108 (1967). 
[7] M. Stone, Phys. Lett. 67B, 186 (1977). 

[8] C. G. Callan and S. Coleman, Phys. Rev. D 16, 1762 (1977). S. Coleman, in: The 
Whys of Subnuclear Physics, edited by A. Zichichi (Plenum, New York, 1979). 

[9] W. H. Miller, J. Chem. Phys. 62, 1899 (1975); W. H. Miller, Adv. Chem. Phys. 25, 
69 (1974). 

[10] E. PoUak, Phys. Rev. A 33, 4744 (1986); E. PoUak, Chem. Phys. Lett. 177, 178 (1986). 

[11] P. Hanggi and W. Hontscha, J. Chem. Phys. 88, 4094 (1988); P. Hanggi and W. 
Hontscha, Ber. Bunsenges. Phys. Chem. 95, 379 (1991). 

[12] H. Grabert, P. Olschowski, and U. Weiss, Phys. Rev. B 36, 1931 (1987). 

[13] M. H. Devoret, D. Esteve, C. Urbina, J. Martinis, A. Cleland, and J. Clarke, in: Quan- 
tum Tunneling in Solids, edited by Yu. Kagan and A. J. Leggett (Elsevier, New York, 
1992). 

[14] R. P. Feynman and F. L. Vernon, Ann. Phys. (N.Y.) 243, 118 (1963). 

[15] H. Grabert, P. Schramm, and G.-L. Ingold, Phys. Rep. 168, 115 (1988). 

[16] U. Weiss, Quantum Dissipative Systems (World Scientific, Singapore, 1993). 

[17] G.-L. Ingold, Ph.D. thesis (Stuttgart, 1988). 

[18] H. Hofmann and G.-L. Ingold, Phys. Lett. B 264, 253 (1991). 

[19] P. UUersma, Physica 32, 27, 56, 74, 90 (1966). 

[20] R. P. Feynman and A. P. Hibbs, Quantum Mechanics and Path Integrals (McGraw-Hill, 
New York, 1965); R. P. Feynman, Statistical Mechanics (Benjamin, New York, 1972). 

[21] L. S. Schulman, Techniques and Applications of Path Integrals (Wiley, New York, 
1981). 

[22] J. Ankerhold and H. Grabert, Physica A188, 568 (1992). 

[23] J. Ankerhold and H. Grabert, Dissipative Quantum Systems with Potential Barrier. II. 

Dynamics Near the Barrier Top, to be published. 
[24] R. F. Grote and J. T. Hynes, J. Chem. Phys. 73, 2715 (1980); P. Hanggi and F. 

Mojtabai, Phys. Rev. A 26, 1168 (1982). 
[25] H. Grabert, U. Weiss, and P. Talkner, Z. Phys. B 55, 87 (1984). 
[26] P. G. Wolynes, Phys. Rev. Lett. 47, 968 (1981). 



28 



FIGURES 

FIG. 1. Diagonal part gfi(0,g) of the form factor of the stationary flux solution (102) as a 
function of the scaled coordinate q for ohmic damping with 7 = 3 and different values of the 
inverse temperature 9. The dashed line corresponds to ^ = 0.1, the dotted-dashed line to ^ = 0.5, 
and the solid line to ^ = 3.0. 
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